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giving  rise  to  the  non-Debye  relaxation. We  suggest  that  the  microscopic 
mechanism  underlying  this  TDTR  is  another  aspect  of  the  heat  bath  with  which 


the  relaxing  system  is  in  contact,  not  considered  before.  The  nature  of  the 
low-energy  level  spaclngs  of  the  heat  bath  determine  the  long  time  relaxation 


process.  Alternative  procedures  of  arriving  at  the  same  TDTR  are  also 
discussed. 
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The  slower  than  Debye  relaxation  In  condensed  matter  Is  well-known  as  is 
evidenced  by  this  topical  meeting  itselfl  In  this  paper,  we  develop  in  detail 
the  theoretical  ideas  proposed  earlier  by  one  of  us*  some  time  ago,  with  the 
hope  that  some  of  the  obscurities  in  that  paper  are  removed  and  our  current 
thinking  on  the  subject  is  explained  in  some  detail.  Ngai  has  given  a 
detailed  account  of  the  experimentally  observed  universal  form  for  the 
relaxation  of  any  physical  quantity,  viz. 

exp  -  a  (t/tp)k,  a>0,  0<b<l  (1) 

with  t  p,  a  characteristic  time  In  the  system,  under  diverse  physical 
conditions.  That  the  form,  Equ.(l)  is  so  universally  found  is  intriguing  and 
we  seek  a  theoretical  explanation  of  It  which  is,  as  far  as  possible,  also 
universal  and  model -Independent.  That  the  relaxation  be  bounded  by  a  function 
of  the  form  (1)  can  be  deduced  by  the  mere  requirement  of  the  semi  boundedness 
of  the  spectrum  of  the  Hamiltonian  describing  the  relaxing  system  and  its 
attendant  environment  and  the  square  integrability  of  the  relaxation  function 
itself,  by  using  the  Paley-Wiener  (PW)  theorem,  was  recently  proposed3.  This 
observation  gives  rise  to  the  Paley-Wiener  criterion  for  relaxation  which  has 
the  strong  Implication  that  Indeed  the  usual  Debye  form  with  b*l  in  Equ.(l)  is 
ruled  out  because  then.  It  would  Imply  the  spectrum  of  the  Hamiltonian  is  not 
bounded  from  below,  which  Is  not  physically  acceptable.  The  aspect  of  square 
Integrability  of  the  relaxation  function  Is  new  and  in  the  absence  of  no 
Information  as  to  the  nature  of  the  long  time  behaviour,  this  requirement  does 
not  seem  unpalatable.  The  quantum  statistical  evolution  of  a  system  Is 
controlled  by  the  Liouville  operator  given  the  initial  state  of  the  system  but 


no  other  statement  made  as  to  the  nature  of  the  time-dependence,  except 
that  for  large  times  the  system  evolves  Into  a  statistical  equilibrium  state. 
No  comment  can  be  made  In  general  as  to  how  this  state  Is  approached  with  the 
present  set  up  of  the  equations  and  conditions  Imposed.  If  one  adds  the 
physical  requirement  of  the  PW,  then  one  can  make  detailed  comments  about  the 
time  dependent  of  approach  to  equlllthlum.  This  PW  criterion  Is  thus  a 
postulate  which  Is,  as  of  our  present  understanding,  a  general  statement. 
Independent  of  the  actual  model.  But  In  actual  fact,  the  bound  Is  attained 
experimentally  and  it  is  this  aspect  that  was  addressed  to  in  Ref.  1.,  where 
It  was  postulated  that  the  long-time  behaviour  (viz.  longer  than,  say,  10"10 
sec)  is  controlled  by  low  energy,  many-body  excitations  of  the  system  with 
energies  less  than,  say,  10+*®  Hz.  This  leads  us  to  the  discussion  of  what  is 
"relaxation"  in  condensed  matter  systems.  One  considers  a  physical  property 
such  as  electric  or  magnetic  or  thermodynamic  or  structural  etc.  and  examines 
how  this  property  behaves  as  time  evolves.  We  think  of  these  properties  as 
being  described  by  the  respective  "species"  such  as  dipoles  for  dielectric 
relaxations,  spins  for  magnetic  problems  etc.  The  evolution  at  long  times 
(  >  10“*®  sec)  is  the  regime  that  Is  being  considered  here.  This  In  turn 
Implies  that  we  are  probing  energy  levels  of  the  systems  of  order  h/10“*°  sec 
or  10“*  °K.  Had  we  considered  shorter  time  scales  such  as  10“13  sec,  the 
corresponding  energy  levels  probed  are  of  order  300°K,  corresponding  to 
excitations  of  the  system  via  phonons  etc.  But  at  energies  such  as  10“*  °K, 
one  has  no  known  excitations  and  one  must  seek  newer  low  energy  mechanisms 
which  drive  the  system  to  relax.  Moreover,  this  process  is  necessarily 
"classical"  Involving  such  low  energies  because  equilibrium  thermal  energies, 

kT,  at  which  experiments  are  conducted  are  almost  always  much  larger  than  10“* 
°K.  We  ask  ourselves  what  these  excitations  could  be  and  how  can  we  describe 


them?  Ordinary  quantum  statistical  mechanical  description  must  accomodate 
these  low  energy  processes  in  a  crucial  way  in  order  to  be  able  to  give  nse 
to  time-dependence  of  relaxation  of  the  form  given  by  Equ.(l).  From  the 
Paley-Wiener  argument,  we  noted  earlier  that  the  pure  Debye  exponential  decay 
is  not  acceptable  because  the  Hamiltonian  of  the  system  will  not  then  be 
bounded  from  below.  This  difficulty  Is  removed  if  we  considered  the  coupling 
of  the  "species"  to  the  low  energy  excitations  of  a  Heat  Bath  in  which  the 
given  system  is  immersed.  The  Heat  Bath  is  assumed  to  be  a  larger  system 
which  is  itself  unaffected  by  the  presence  of  the  system  under  consideration 
but  causes  the  system  to  relax  to  equilibrium  which  is  characterized  by  the 
temperature  of  the  heat  bath.  The  interaction  between  the  system  and  the  bath 
is  such  that  the  relaxation  comes  about  without  anything  happening  to  the  bath 
Itself.  Conventional  Hamiltonian  description  sets  a  basic  energy  scale  and  in 
the  case  of  a  heat  bath,  it  can  only  be  a  continuous  energy  spectrum  with  no 
apparent  discrete  energy  scale  because  the  bath  is  considered  very  large 
compared  to  the  system  with  which  It  is  in  contact  and  only  the  average  of  the 
bath  Hamiltonian  is  needed  to  specify  the  temperature  to  which  the  system 
ultimately  goes  to.  It  therefore  appears  abundantly  clear  that  the  bath 
Hamiltonian  cannot  be  specified  in  detail  but  yet  what  is  needed  in 
determining  the  long  time  relaxation  phenomenon  Is  the  small  energy-level 
spacings  (of  order  of  10"4  ev  or  lower)  of  the  spectrum  of  the  heat  bath. 

We  have  thus  motivated  the  use  of  a  non-specific  large  system  whose 
average  spectral  properties  are  only  needed  in  developing  a  theory  of 
relaxation  phenomenon  in  condensed  matter  systems.  We  note  at  this  stage,  the 
description  of  the  heat  bath  by  the  Wigner  Random  Matrix  Hamiltonian  theory 
seems  admirably  suited  with  our  requirements  stated  above.  We  now  proceed  to 
outline  our  theory  of  relaxation.  This  involves  three  steps:  (1)  derivation 


of  the  relaxation  function;  (2)  Introduction  of  the  jump-rate  of  transition; 
and  (3)  the  calculation  of  the  jump-rate. 

We  consider  some  elementary  entitles  as  being  the  actors  In  the  relaxation 
process.  They  flip  between  equivalent  energy  states  (energy  degenerate)  and 
are  thought  to  be  Immersed  in  a  bath  system  described  by  the  Wigner  Random 
Matrix  Hamiltonian.  We  use  the  Gaussian  Orthogonal  Random  Matrix  Hamiltonian 
(GOE)  because  we  consider  the  heat  bath  to  be  time-reversal  Invariant  and  also 
that  it  has  a  definite  average  energy  so  that  one  may  employ  a  canonical 
ensemble  picture  to  describe  It  with  a  fixed  temperature  associated  with  it. 

In  equilibrium  then,  the  actors  acquire  this  same  temperature,  T. 


Step  1:  Derivation  of  the  Relaxation  Function 
Let  1,2  represent  (for  simplicity  we  consider  a  two-equivalent -state 
system)  the  two  equivalent  system  states  and  let  pj(t)  and  p2(t)  be  their 
occupancies  at  an  instant  t.  The  Interaction  of  these  system  states  with  the 
heat  bath  Induces  transitions  between  these  two  states.  We  are  Interested  in 
computing  the  relaxation  of  this  system  which  is  merely  the  correlation  of 
P(t)  =Pi(t)  -  P2 (^ )  at  two  different  times.  We  are  thus  Interested  in 

*(t)  =  -  «P(tQ)  P(t0-t)»  (2) 

« — »  Is  the  quantum  statistical  thermal  average  over  the  entire  bath  plus 
system  hamlltonlan4.  Here  t0  Is  some  Initial  time.  We  are  Interested  In  the 
long  time  (t  +  •)  behaviour  of  this  function.  Thus,  for  frequencies 
w«10  Hz,  which  Is  equivalent  to  a  temperature  of  1/2  °K  and  below,  we  may 
make  the  approximation  kT»^wand  so  only  the  classical  limit  of  Equ.(2) 
suffices.  In  this  limit,  Equ.(l)  becomes 


(3) 


♦(t)  -  -  8<P  (tQ)  P  (t0-t)>Q 

where  < — >0  Is  the  classical  average  over  the  equilibrium  distribution 
function  and  g«l/kT.  P  here  denotes  time  derivative.  It  is  important  to 
stress  here  that  in  the  long  time  limit,  the  classical  limit  of  (2)  suffices 
because  kT>>^u*»  of  the  system. 


Step  2:  Introduction  of  the  Jump-transition  Rate 

In  the  classical  limit  then,  it  suffices  to  consider  a  rate  equation  for 
the  populations  for  Pj(t)  and  p2(t)  and  deduce  therefore  a  form  for  the 
relaxation  function,  (3).  Defining  W(t)  to  be  the  general  time -dependent  jump 
transition  rate  for  flipping  from  state  1  to  2  and  vice  versa,  we  have,  the 
rate  equations 


0x(t)  =  -  W(t)  (pj (t )  -  p2(t )) 


02(t)  *  -  W(t)  (p2(t)  -  Pl(t)) 


(4) 


where  we  have  taken  into  account  the  fact  that  the  flipping  occurs  between  the 
two  equivalent  states  only  and  that  there  is  no  loss  of  species  in  the  process 
i.e.,  “0.  Thus,  we  find  at  once 


P(t)  -  -  2  W(t)  P(t) 


(5) 


where  P(t)  ■  Pj(t)  -  p2(t)  introduced  before  in  Step  1.  Thus,  we  obtain, 
P(t)  •  P0  expf-z/*  H(f)dt')  (6) 


where  t  P0  are  the  values  that  the  system  takes  as  time  progresses, 
the  relaxation  function,  (3),  has  the  form 

*(t)  =  2  e  P02  W(t)  exp  (-2 £  W(t')dt') 


And  so 

(7) 


It  should  be  observed  that  If  W(t)  Is  taken  to  be  a  constant.  Independent 
of  time,  say  W0,  then  we  obtain  the  classical  Debye  result: 

+(t)  =  2  0  P  2W  exp (-2  W  t)  (8) 

Debye  00  0 

and  the  relaxation  time  Is  t  «  1/2Wq.  In  obtaining  the  rate  (or  master) 
equation,  (4),  one  often  goes  through  a  detailed  analysis  of  the  density 
matrix  with  Markofflan  assumptions  concerning  the  process  that  are  going  on  In 
the  system  and  the  constant  WQ  Is  the  time-independent  transition  rate  (TITR) 
given  by  the  Golden  Rule  of  Quantum  Mechanics.5  It  should  therefore  be 
stressed  that  the  steps  I  and  II  above  are  really  stemming  from  more 
fundamental  considerations  even  though  we  have  presented  them  here  In  a 
simpler,  physical  way.  The  question  then  Is  to  consider  the  evaluation  of 
W(t)  In  greater  detail  than  has  been  done  before. 

Step  III:  Calculation  of  W(t) 


The  mechanism  that  gives  rise  to  W(t)  Is  the  Interaction  of  the 
equivalent  states  Introduced  In  Step  I  with  low  lying  manybody  states  with 
energies  less  than  say,  lO^Hz,  In  order  to  determine  the  long  time  decay 


phenomena.  In  doing  this  computation,  the  spontaneous  emission  and  absorption 
processes  are  both  taken  Into  account.  In  fact,  this  calculation  Is  very 
similar  to  the  phonon  side  band  problem  except  that  In  the  present  context  we 


have  to  work  with  the  Wlgner  system  Instead  of  the  phonons.  The  result  Is 

W(t)  =  WQ  exp(-$<t))  (9) 

where  ♦(t)  Is  the  time  response  of  the  heat  bath  (or  Wlgner  system),  taking 
care  to  remember  that  the  energies  Involved  here,  i\  w,  are  all  much  less  than 
kT,  the  thermal  energy,  so  that  the  associated  probabilities  for  both  emission 
and  absorption  are  just  constants,  a  factor  of  1/2.  Moreover,  here  4>(t)  Is 
real  with  no  threshold  effects. 

In  order  to  derive  Equ.(9)  and  an  explicit  form  for  $(t)  we  proceed  as 
follows.  We  consider  a  relaxing  system  l.e.,  after  a  certain  relaxation 
time  t,  the  system  Is  thought  to  be  In  a  relaxing  state.  Given  such  a  t,  we 
define  destruction  and  creation  operators  R  and  R*  associated  with  relaxation 
such  that 

R|>  *  0  In  the  relaxing  state  (10) 

where  |>  refers  to  a  pseudo  "vacuum"  of  the  relaxing  state  and  the  heat  bath 
with  which  the  system  Is  in  contact.  Then  R*(o)|>  is  the  given  relaxing  state 
at  an  Initial  reference  time  (this  Is  the  time  the  system  had  begun  to  relax!) 
which  we  denote  here  as  the  zero  of  our  reckoning  of  time.  We  then  ask  how 
this  state  evolves  In  time.  This  time  evolution  Is  governed  as  usual  by 

R*  (t ) | >  *  e*Ht  R*(o)e"^Ht|>  (11) 

where  H  Is  the  Hamiltonian  describing  the  relaxing  system,  the  heat  bath  with 
which  It  is  In  contact,  and  the  Interaction  between  them.  This  we  denote  by 
H0  +  V.  We  are  thus  led  to  the  calculation  of  the  quantity 


b(t)  =  <| R (t )  Rf  (0) | > 


=  <|eiHt  R(o)  e"iHt  Rf  (o)|>  (12) 

Since  R*  (o)|>  represents  the  given  relaxing  state, 
e“iHt  Rt(o)|>  =  Rf(o)  *xp-i(H0+V)Rt 

where  (HQ+  V)R  stands  for  the  part  of  the  Hamiltonian  associated  with  the 
relaxing  state  under  study.  Since  we  are  here  for  simplicity,  considering 

only  one  relaxing  state,  R(o)  R*(o)  =  1,  and  moreover  <|e+i^o+v^t=  ^e^o* 
because  <|  by  virtue  of  its  definition  contains  no  interaction  between  the 
relaxing  state  and  the  heat  bath.  Thus, 

b(t)  =  <|exp  (i  H0t)  (exp  -i(HQ+V)Rt)  l>  (13) 

The  state  R  in  the  above  does  not  contain  any  index  because  we  will  invoke  the 
"randomness"  of  the  Heat  Bath  Hamiltonian  to  eliminate  reference  to  any 
special  state.  With  this  understanding,  we  may  drop  the  subscript  R  in 
Equ.(13)  and  compute 

U(t)  •  e'Ve-'IV'')'  <14> 

by  using  standard  cluster  expansion  techniques*  The  result  Is 

b(t)  »  exp  {♦(t )}  (15) 


where 


(15a) 


We  have  thus  reduced  the  problem  to  consideration  of  the  bath  states 
only.  Remembering  that  the  bath  states  are  treated  in  a  canonical  ensemble 
framework,  and  since  $.(t)  involves  expectation  value  of  the  interaction 
between  the  system  and  the  bath,  it  may  be  dropped  from  further  consideration 
and  consider  only  <t>2(t).  The  other  higher  order  terms  may  be  dropped,  partly 
because  of  the  weak  coupling  between  the  bath  and  the  system  and  partly 
because  it  would  involve  the  multiple  energy  level  distributions  which  may  be 
expected  to  be  of  smaller  significance  compared  to  the  first  nontrivial 
contribution,  <t>2  .  Thus 


V  ,  2 

♦2(t)  -  -V2Z  y-r^r-To  h  Q-e  n  (Ei  "Ej  h  +  *5  (1-e' 
i,j  Ui'\j^ 
i*j 


The  factors  V2  in  {....}  appear  from  equal  occupancies  of  the  state  of 

the  bath.  The  sum  over  i,j  can  be  written  as 


♦2(t> .  -  j 


l»..i  l2>  < 


1-coset. 


(17) 


We  now  argue  that 


V2 


E.-E 


| V i j  j 2 )  is  Independent  of  e  because  | j | 2  is 
e 


a  random  quantity  uncorrelated  with  e  except  perhaps  when  there  are  strong 
correlations  which  may  occur  when  the  energies  are  high  such  that  the  modes 
are  well  localized.  This  aspect  Is  the  basic  strength  of  the  approach  of  the 
model  because  the  "randomness"  now  has  been  separated  out  of  the  main  problem 
so  that  a  universal  result  Is  obtained  without  the  knowledge  of  the  random 
distribution  of  |V^j|^w1th  E^-Ej=  e  •  Thus,  Introducing 

ec 

£  fol  Vij  2)+f  de  |V|2  N(e) 

1*j  0 

E.-E.  =  e 
i  J 

we  obtain, 

♦o(t)  =  /c  de  N(e)  | v| ^{l^COSet }  (18) 

Jo  e 

The  cut  off  ec(»10^Hz)  Is  Introduced  here  such  that 

| Vj2  N  (e)  =  en,  (0  <  n  <  1)  (19) 

Is  violated  beyond  e  *  ec#  Here  we  have  invoked  the  Wigner  result  concerning 
the  level -spacing  distribution  of  a  GOE.  This  final  result  then  leads  to  the 
universal  form*  of  relaxation  given  In  Equ.(l)  with  b*l-n. 

The  above  arguments  are  at  a  physical  level,  lacking  perhaps  the  rigor  of 
mathematics.  Our  aim  was  to  give  a  simple  physical  derivation  of  Equ.(l). 
There  Is  a  completely  heuristic,  mathematical  way  of  arriving  at  such  a  result 
also  by  Invoking  the  mathematical  Idea  of  "fractional  Brownian  Notion"  (and 
the  fractal  times)  of  Mandelbrot  and  van  Ness^.  This  procedure  would  Imply 


that  relaxation  process  may  be  thought  of  as  basically  a  generalized  diffusion 
process.  The  conventional  diffusion  process  using  a  Langevin  equation  with 
white  noise  stochastic  process  leads  to  the  Debye  exponential  relaxation.  If 
one  believes  in  the  diffusion  mechanism  as  underlying  relaxation  processes, 
then  one  is  forced  to  modify  the  white  noise  stochastic  picture  of  diffusion 
by  the  "fractional  Brownian  noise"  so  that  one  obtains  a  non-Debye  relaxation 
of  the  form  given  by  Equ.(l).  But  we  know^  that  not  all  physical  relaxation 
processes  can  be  related  to  Brownian  motion  and  so  this  model  suffers  from  a 
lack  of  physical  basis  in  spite  of  the  fact  that  one  may  be  able  derive  non¬ 
exponential  decay  of  the  form  given  by  Equ.(l).  This  is  an  interesting 
mathematical  aspect  of  the  non-Debye  relaxation.  A  purist  may  object  to  this 
approach  to  relaxation  because  the  basic  premise  here  is  that  the  underlying 
rate  equation  is  a  Langevin  type  equation  with  fractional  brownian  noise,  to 
which  one  has  no  answer  except  to  say  that  it  seems  to  lead  us  to  the  observed 
result.  There  are  attendant  changes  in  the  Einstein  relation  etc.  that  ensue 
such  a  framework  and  the  reader  may  refer  to  Ref.  8  for  a  discussion  of  it. 

There  is  another  level  at  which  the  present  problem  may  be  attacked. 

This  is  perhaps  the  most  complete  theory  but  even  here  certain  characteristics 
of  the  heat  bath  have  to  be  introduced  at  a  particular  point.  We  present  this 
approach  for  the  complete  non-believer  with  the  hope  that  this  approach  may  be 
appealing  to  both  physically  and  mathematically  minded  person.  In  this 
approach  one  begins  with  a  density  matrix  for  a  system,  heat  bath,  and  their 
mutual  interaction  with  the  explicit  knowledge  that  the  heat  bath  Is  much 
larger  than  the  system  In  contact  with  It  and  their  mutual  interaction  Is  weak 
such  that  the  leading  order  calculation  suffices.  Moreover  one  constructs  the 
system  density  matrix  by  taking  a  trace  over  the  bath  so  that  the  result  Is 
free  of  all  bath  variables.  The  difference  between  this  approach  and  a 


similar  approach  found  In  the  literature®  In  the  derivation  of  the  master 
equation  for  the  system  density  matrix  Is  that  In  our  procedure,  while  we  make 
the  Markoff  approximation,  we  keep  the  time-dependence  of  the  transition 
rate.  By  so  doing,  we  bring  In  details  of  the  bath  Hamiltonian.  Also, 
Invoking  the  Wlgner's  random  matrix  Hamiltonian  as  a  description  of  the  heat 
bath,  which  seems  most  appropriate  for  describing  the  unknown,  large  system, 
we  see  that  this  program  makes  the  derivation  of  an  equation  such  as  Equ.(5) 
including  an  expression  for  W(t)  quite  acceptable,  on  a  formal  level. 

In  this  case,  we  do  not  invoke  the  PW  theorem  but  we  are  in  comformity 
with  It.  We  use  the  known  Information  about  the  system  plus  heat  bath  at  the 
initial  time  where  It  Is  given  that  the  bath  is  in  a  canonical  ensemble  and 
the  system  Is  In  some  preassigned  system  density  matrix  and  the  above  program 
answers  the  question  of  ultimate  time  evolution  of  the  system  to  an 
equilibrium  ensemble  state.  In  so  doing,  unlike  In  the  equilibrium 
statistical  mechanical  theory  where  the  system  relaxes  to  an  equilibrium  state 
with  the  same  temperature  as  the  heat  bath,  which  Is  merely  the  average  of  the 
heat  bath  hamlltonlan,  the  long  time  relaxation  behaviour  depends  on  the 
nature  of  the  bath  spectrum  for  small  energy  level  spaclngs.  Using  the 
results  of  the  Gaussian  Orthogonal  ensemble  theory  of  random  matrix 
Hamiltonian  as  was  done  In  Ref.  1,  we  deduce  the  universal  relaxation  of  the 
form  given  of  Equ.(l).  For  a  discussion  of  this  theory  but  In  the  context  of 
quantum  mechanics  of  the  system  plus  bath  system,  one  may  refer  to  the  paper 
of  Rajagopal  and  Wlegel.9 

In  the  derivation  of  the  universal  relaxation  function,  Eq.(l)  given 
earlier,  the  bath  states  take  a  subservient  role.  The  relaxation  species  Is 
the  center  of  attention.  Its  relaxation  rate  Is  calculated  and  the  bath 
states  come  Into  the  picture  only  In  how  th^y  modify  Its  relaxation  rate. 


Stimulated  by  an  Idea  of  Palmer  and  Anderson10,  we  shall  give  here  another 
derivation  of  the  universal  law  (1)  while  treating  both  the  relaxation  species 
and  the  bath  states  altogether  on  an  equal  footing.  That  Is,  we  start  by 
writing  down  the  relaxation  rate  as11 

W(t)  =  W0e-(£  -  T  AS)/kT  (2°) 

Here  WQ  exp(-E/kT)  represents  the  transition  rate  of  the  relaxation  species  If 
Its  Interaction  with  the  bath  states  Is  Ignored.  Thus  E  Is  an  activation 
energy  if  Its  value  is  nonzero.  WQ  is  the  "attempt  frequency"  of  the  primary 
species.  The  quantity  AS  is  the  entropy  of  the  bath  states  measured  from  its 
maximum  value.  Initially  before  the  commencement  of  relaxation,  the  bath 
states  are  in  equilbrium  and  Its  entropy  must  have  a  maximum  value,  Smax. 

After  relaxation  has  commenced  at  t=0,  the  interaction  between  the 
relaxing  species  and  the  bath  states  causes  excitation  and  deexcitation  of  the 
bath.  We  calculate  the  probability  that  the  bath  remains  in  the  same  state 
after  the  Interaction  V  Is  turned  on.  In  the  lowest  order  perturbation 
theory.  This  probability  Is 
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because  we  assume  that  the  bath  states  are  perturbed  very  little  by  the 
Interaction  with  the  system.  If  the  Interaction  V  were  neglected,  p^  *  1  as 
It  should,  by  construction.  Then,  calculating  the  entropy  change  due  to  this, 
we  have,  to  leading  order  In  the  Interaction, 


(22) 


In  the  above,  1 , j  refer  to  a  group  of  states  of  the  bath  which  form  the  bath 
ensemble. 

Following  the  same  arguments  as  before,  (see  Equ.  17,18,19)  we  find 


1  E  AS 

between  the  observed  xp  and  the  microscopic  tq  =  (Wq  exP(pf  -  •=£>)). 

Obviously  the  present  approach  and  the  earlier  one  are  closely  related. 
Both  have  the  Infrared  divergent  bath  states  excitation  at  their  root.  The 
difference  between  the  earlier  approach  and  the  present  one  Is  that  the  entire 
system,  relaxation  species  plus  the  bath  states  are  considered  at  the  same 
time.  It  Is  not  surprising  that  we  obtain  the  same  results  Independent  of  the 
approach.  However  the  present  approach  which  lets  Its  bath  states  play  a  more 
active  role  enables  further  Interpretation  of  the  physical  processes  that 
occur  In  relaxation  phenomena  In  condensed  matter.  Relaxation  of  the  species 
will  lead  to  an  Increase  In  entropy  with  time  as  all  Irreversible  processes  do 
according  to  the  second  law  of  thermodynamics.  What  Is  Interesting  here  Is 
that  not  only  the  relaxation  Is  non-exponential  In  time  and  have  the  universal 
form  of  Eq.(l),  but  also  this  relaxation  Is  accompanied  by  a  decrease  In 
entropy  of  the  bath  states. 

In  closing  we  would  like  to  discuss  another  derivation  of  the  universal 
relaxation  function  which  Is  totally  different  In  spirit  from  the  mechanistic 
models  given  earlier.  Given  the  universal  form  of  the  macroscopic  relaxation 
function  In  Eq.(l),  one  Is  led  to  consider  a  microscopic  mechanism  underlying 
relaxation  phenomenon  In  general.  We  may  think  of  the  relaxing  body  as  being 
composed  of  some  species  of  relaxing  entitles  such  as  electric  dipoles  In  a 
dielectric,  charge  carriers  In  a  semiconductor,  etc.,  each  of  which  had 
Identically  distributed  energy  variables  with  the  same  energy  distribution. 

The  macroscopic  energy  distribution  function  p(e)  associated  with  the 
relaxing  body  for  low  energies  Is  the  limit  distribution  of  normalized  sums  of 
the  microscopic  energy  variables.  A  similar  procedure  Is  Involved  In 
statistical  mechanics  for  deriving  equilibrium  ensemble  distributions.  The 
resulting  macroscopic  density  distribution  p(e)  necessarily  obeys  the 


relationship 
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p(ae+b)  =  (a^/a)  f  de'pfa^e-e'J+b^pCa^'+bg)  (2) 

•■CD 

for  every  set  of  parameters  a^  >  0,  b-,  with  the  corresponding  a  >  0,  b  related 
to  the  parameters  a^,  b^.  The  relaxation  function  Is  proportional  to  the 
modulus  square  of  the  Fourier  transform,  jc(t)|^,  of  the  energy  distribution 
function  p(e).  For  distributions  obeying  Eq.‘(2)  this  is  necessarily  of  the 
form^ 


| c(t ) | 2  =  e  a>0,  0<b<2  (3) 

If  In  addition,  the  spectrum  of  the  Hamiltonian  is  semlbounded,  then  the 
Paley-Wiener  theorem  applies  and  the  values  of  b  would  have  to  be  in  the 
smaller  range  0<b<l,  leading  us  again  to  Eq.(l).  For  such  microscopic  models, 
we  see  that  the  Paley-Wiener  bound  becomes  exact. 

We  have  thus  given  here  a  variety  of  general  arguments  at  different 
levels  to  obtain  the  non-Debye  relaxation  in  condensed  matter  properties. 
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before  publication. 
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